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■ ABSTRACT 

OO ' 

I ■ We develop a model to implement metal enrichment in a cosmological 

context based on the hydrodynamical AP3MSPH code described by Tissera, 
C*") ■ Lambas and Abadi (1997). The star formation model is based on the Schmidt 

law and has been modified in order to describe the transformation of gas into 
stars in more detail. The enrichment of the interstellar medium due to super- 
f—^ novae I and II explosions is taken into account by assuming a Salpeter Initial 

t~^. ■ Mass Function and different nucleosynthesis models. The different chemical 

elements are mixed within the gaseous medium according to the Smooth Par- 

■ ticle Hydrodynamics technique. Gas particles can be enriched by different 

neighbouring particles at the same time. We present tests of the code that 

& ■ assess the effects of resolution and model parameters on the results. We show 

that the main effect of low numerical resolution is to produce a more effec- 
Q ■ tive mixing of elements, resulting in abundance relations with less dispersion. 

| We have performed cosmological simulations in a standard Cold Dark Mat- 

ter scenario and we present results of the analysis of the star formation and 
\ chemical properties of the interstellar medium and stellar population of the 



^1- 



simulated galactic objects. We show that these systems reproduce abundance 
ratios for primary and secondary elements of the interstellar medium, and the 
jJJ ■ correlation between the (O/H) abundance and the gas fraction of galaxies. We 

find that star formation efficiency, the relative rate of supernovae II to super- 
novae I and life-time of binary systems as well as the stellar nucleosynthesis 
model adopted affect the chemical properties of baryons. We have compared 
the results of the simulations with an implementation of the one-zone Simple 
Model, finding significant differences in the global metallicities of the stars 
and gas as well as their correlations with dynamical parameters of the sys- 
tems. The numerical simulations performed provide a detailed description of 
the chemical properties of galactic objects formed in hierarchical clustering 
scenarios and proved to be useful tools to deepen our understanding of galaxy 
formation and evolution. 

Key words: cosmology: theory - galaxies: formation - galaxies: evolution - 
galaxies: abundances. 

1 INTRODUCTION 

In recent years, our knowledge of the high redshift Universe has increased dramatically 
allowing us to start constructing a picture of how the different morphological types evolve 
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(e.g. Steidel and Hamilton 1992; Madau 1995; Lilly et al. 1995; Cowie et al. 1996; Ellis et al. 
1996; Steidel et al. 1998). In particular, the comoving star formation history first depicted by 
Madau et al. (1996) has resulted in a very useful way of quantifying galaxy evolution, under 
certain hypothesis. Since star formation directly translates into metal enrichment of the 
interstellar medium (ISM), observations of these two processes help astronomers to constrain 
models of the formation of the structure. In particular, Lyman alpha forests, damped Lyman 
alpha systems and Lyman alpha break galaxies have all contributed to estimate the cosmic 
metal ejection and star formation rates (e.g., Lu et al. 1996; Pettini et al. 1997). However, 
we are still far from drawing a consistent picture since several points remain to be clarified. 
In this respect, the integration of the dust corrected cosmic star formation history of the 
Universe up to z = can account for the entire stellar mass content in spirals and spheroids 
at present times. However, when we look at the metal content at z ~ 2.5, only 10% of 
what is expected is actually measured (Renzini 1998). Regarding the chemical properties of 
galaxies, most of the available observations are restricted to the Galaxy, and in particular, 
to the solar neighbourhood (e.g., Edvardsson et al. 1993; Gratton et al. 1996; Rocha-Pinto 
et al. 2000). Extragalactic observations of HII regions provide information on the chemical 
content in other galaxies (e.g., Pagel 1992; Garnett et al. 1995; Kennicutt and Garnett 1996; 
Kobulnicky and Skillman 1996). 

Chemical evolution models are a useful tool to attempt to study the physical processes 
that might determine the chemical characteristics of the different galaxies. Detailed mod- 
elling can be found in studies of the Galaxy (e.g., Ferrini et al. 1992; Tossi 1996; Chiappini, 
Matteucci and Gratton 1997 and references therein). These analytical models carefully treat 
stellar evolution. However, they cannot account for dynamical evolution, and certainly, not 
for hierarchical clustering. From a numerical point of view, metallicity enrichment mech- 
anisms have been implemented in hydrodynamical simulations in different ways. The first 
attempts were done by Larson (1975, 1976). Other implementations came after using diverse 
techniques such as chemodynamical models (e.g., Burkert et al. 1992; Samland et al. 1997) 
which describe in more detail the interstellar medium evolution. On the other hand, this 
approach uses a very simple prescription for galaxy formation where the dark matter halo 
is either not considered at all or assumed to be static. However, there are strong evidences 
from observations of normal spirals that the dark matter is dynamically important within 
the luminous radius (e.g., Bottema 1992; Courteau, de Jong & Broeils 1996; Rhee 1996), 
and that the response of the dark matter to the presence of baryons affect their evolution 
(Blumenthal et al. 1996; Tissera & Dommguez-Tenreiro 1998) and, as a consequence, the 
star formation process (Mihos & Hernquist 1996; Barnes & Hernquist 1996; Tissera 2000; 
Navarro & Steinmetz 2000). A first approach to include chemical evolution in Smooth Parti- 
cle Hydrodynamical (SPH) code is described by Steinmetz and Miiller (1994, 1995) followed 
by Raiteri, Villata and Navarro (1996). These works run prepared-cosmological initial condi- 
tions where the formation and evolution of only one object was studied. These models have 
shed light on some physical mechanisms that may control metallicity in galaxies, and have 
also shown that hierarchical clustering scenarios may form galactic objects that resemble 
Milky Way-type galaxies from a chemical point of view. However, the SPH models with 
chemical implementations hitherto published are restricted to the Galaxy. It would be very 
important for the study of the formation and evolution of galaxies to be able to analyse 
a variety of galactic objects with different evolutionary histories. Hydrodynamical cosmo- 
logical simulations are more suited to tackle this problem since the non-linear evolution of 
the matter is naturally accounted for, and, physical processes can be more consistently im- 
plemented. These models can provide coherent well-described environments for all objects 
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and a complete record of their formation and evolution. The drawback is the treatment of 
numerical resolution effects. 

In this paper, we concentrate on the description of the chemical model that has been 
implemented in a cosmological hydrodynamical AP3M. First results on galaxy formation, 
global chemical properties and comparison with observations are reported. The detailed 
analysis of the chemical properties of the ISM and stellar population in galaxy-like objects 
such as abundance gradients, age-metallicity relations, etc., are given by Tissera et al. (2000, 
hereafter Paper II). 

This paper is organized as follows. Section 2 describes the star formation process and the 
metal enrichment implementation. In Section 3 we assess the performance of the chemical 
model. In Section 4 we present the first results on galaxy formation and chemical evolution. 
Section 5 summarizes the results. 



2 NUMERICAL MODEL 

We use a cosmological numerical code based on the SPH technique described by Tissera et 
al. (1997). The SPH algorithm has been coupled to the AP3M gravitational code (Thomas 
and Couchman 1992), in order to follow the gravitational and hydrodynamical evolution of 
particles within a cosmological context. For the sake of simplicity, and as first step in the 
development of this chemical cosmological code, simulations have been run considering that 
the gas cools down by using the approximation given by Dalgarno and McCray (1972). 

2.1 Star formation 

In general, the modelling of star formation (SF) in SPH simulations is based on the Schmidt 
law and a series of hypothesis to select suitable SF regions. However, how it is actually 
implemented within each particular code depends on the different authors (e.g., Katz 1992; 
Navarro and White 1994; Tissera et al. 1997). In this paper, we describe a modified version of 
the SF implemented by Tissera et al. (1997) in order to be able to track the transformation 
of gas into stars in more detail, within a given gas particle. 

We include the SF algorithm as follows. Gas particles are eligible to form stars if they 
are cold (T < T*) and satisfy a density criterium: p gas > p crit . This density criterium arises 
by requiring the cooling time of a gas particle to be smaller than its dynamical time. The 
critical temperature is taken as the minimum provided value by the cooling functions 
(T* ~ 10 4 K). Finally, gas particles have to be part of a collapsing region. This requirement 
is imposed by selecting particles in a convergent flow (x/.v < 0). When a gas particle satisfies 
all these conditions, star formation occurs according to the Schmidt law, 

dpstai Pgas /, \ 

^r = c i7' [) 

where c is the star formation efficiency and is a characteristic time-scale assumed to be 
proportional to the dynamical time of the particle (t* = t^ yn = (3n /lQGp gas ) 1 ^ 2 ). Then, each 
stellar mass formed in a particle at a certain SF episode is given by A sta r = C At, where 
At is the time step of integration (At = 1.3 x 10 7 yr) and C a new SF efficiency Note that, 
according to equation (1), in order to go from density to mass, the volume of the new-born 
stars has to be assumed. Taking into account that the process we are modelling happens 
within a gas particle and, that according to observations, SF occurs normally in clusters, 
not in isolation, we assume that the volume occupied by the fraction of new-born stars is 
constant for all them. Thus, it can be absorbed in the new constant, C. 



© 2000 RAS, MNRAS 000, [l]-?? 



4 M. B. Mosconi, P. B. Tissera, D. G. Lambas and S. A. Cora 

According to this SF scheme, in a given baryonic particle there could be several SF 
episodes that had occurred at different times. The gaseous mass of the baryonic particle is 
reduced by A sta r until its gas reservoir is depleted. A minimum gas mass equal to 5 % of 
the initial gas mass component is left over in each particle. When a particle reaches this 
minimum mass, it is transformed completely into a star particle and hereafter behaves as 
collisionless matter. Each A star formed can be followed up in time. The number of stars of a 
given mass within A star is estimated by assuming an Initial Mass Function (IMF). Hence, a 
baryonic particle may be formed by a gaseous and stellar components in different proportions 
according to its history of evolution (hereafter, hybrid particle). The stellar component can 
be made up of different stellar populations with different ages and chemical properties, which 
have been formed at different SF episodes. The chemical properties of each stellar population 
(A star ) reflect the chemical state of the ISM at the time of its formation. 

The decoupling between gas and stars depends on the actual rate at which stars are 
formed in each baryonic particle, that, on its turn, depends on each evolutionary path. 
Hybrid particles may introduce some numerical artifacts since the stellar populations follow 
the gas evolution when they should probably not. However, this is not a trivial problem that 
can be easily resolved. An improvement of the decoupling process is under work and will be 
presented in a separate paper. 

2.2 Metal Production 

Metals are produced and ejected to the interstellar medium at the end of the life of stars. 
Most of chemical elements are ejected by Type II supernovae (SNell) except for the iron 
that is mainly produced by Type I supernovae (SNel). Particles are assumed to be initially 
formed by Hydrogen and Helium in primordial abundances (H = 0.75, He = 0.25). The first 
generations of stars with primordial abundances immediately enrich the ISM from where 
the new generations are born. It has to be stressed that we are not including the effects 
of (thermal or kinetic) energy injection into the ISM due to SN explosions. For the sake of 
simplicity we have splitted the treatment of the feedback problem into two stages. Firstly, 
we include chemical enrichment (this paper), and in a second step, we will intend to develop 
an energy feedback model. There have been several attempts to implement energy feedback 
in SPH codes by either injecting thermal (e.g., Katz 1992) or kinetic (e.g., Navarro & White 
1994; Metzler & Evrard 1994; Navarro & Steinmetz 2000; Springel 2000) energies to the gas 
component due to SN explosions. However, these methods are still controversial leaving this 
problem as an open question for galaxy formation. 

Let us now describe the nucleosynthesis prescriptions adopted in this model: 
Type II SNe: We assume that stars more massive than 8 M & end their lives as Type II SNe. 
In order to estimate their number, we adopt a Salpeter IMF with a lower and upper mass 
cut-off of 0.1 and 120 M , respectively. The IMF adopted gives the total number of SNell 
formed in a certain range of stellar masses in a given A star at a certain time. Supernovae are 
thought to eject their whole metal production within a few 10 7 yr. In particular, we assume 
that their life-time is equal to the integration time-step of the simulations. 

For comparison, we resort to both Portinari, Chiosi and Bressan (1998, P98) and Woosley 
and Weaver (1995, WW95) metal ejecta models and follow the evolution of different elements 
according to the information provided by the authors. P98 give metal yields for stars up to 
120 M & , while WW95 assume that stars larger than 40 M & end up their lives as black holes. 
We consider the following elements according to each author: H, He 4 , C 12 , O 16 , Mg 24 , Si 28 , Fe 56 , 
for both P98 and WW95, and also, N 14 , Ne 20 , S 32 , Ca 40 , Zn 62 , for WW95. 

Type I SNe: Following Matteucci and Francois (1989), we assume that supernovae la 
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originate from carbon deflagration in C-0 white dwarfs in binary systems. It is assumed 
that the masses of these binary systems are likely to be in the range 3 — 16 M . The 
adopted nucleosynthesis prescriptions are taken from Thielemann, Nomoto and Hashimoto 
(1993). Type lb SNe are assumed to be half the total number of Type I SNe and to produce 
only iron (fh 0.3 M per explosion). 

The main interest in including SNI events in the models is that they contribute with a 
substantial amount of iron. This element is fundamental to attempt to reproduce several 
observational results (see also Chiappini et al. 1997). The relative rate of SNI to SNII explo- 
sions has been observationally estimated in the solar neighbourhood and from extragalactic 
sources (van der Bergh 1991). These estimations suggest a range of possible values for the 
total SNI rate respect to that of SNII of 2 < SNRII/SNRI < 3.11. 

Binary systems evolve for a certain period of time (tsNi) during which mass is transfered 
from the secondary to the primary star until the Chandrasekhar mass is exceeded and a 
explosion is triggered. It is generally assumed that a fair fraction of SNel will explode after a 
period of tsNi ~ 10 8 — 10 9 yr (Greggio 1996). It is straightforward to estimate when SNel will 
explode and eject metals to the neighbouring region, since the formation time of each A star , 
in a given baryonic particle, is known. The effects of varying both the SNRII/SNRI relative 
rate (hereafter, 9sn) and tsNi are analysed and results are confronted with observations in 
the following Sections. 

To sum up, the free parameters of our chemical model are: star formation efficiency 
(C), the initial mass function (IMF), its upper and lower mass cut-offs, the relative rate 
of different types of supernovae (Osn), the life-time of binary systems that end as SNel 
(tsNi) and the yields. Regarding the IMF, we are going to adopt a constant and unique one 
throughout this paper. 

2.3 Model for metal ejection 

How the chemical elements are distributed and mixed in the interstellar medium is a complex 
problem which is far from being solved from a theoretical point of view (e.g., Tenorio-Tagle 
et al. 1999). Furthermore, numerical simulations do not reach enough resolution to properly 
treat this small-scale mixing process together with the process of galaxy assembly. Hence, 
the ejection and mixing of elements have to be modelled in heuristic way. 

We use the SPH technique to distribute the metals within the neighbouring sphere of 
the j particle where a A star is formed. According to this technique, the value of a certain 
parameter at the location of j particle can be estimated by using the information storaged in 
its neighbourhood, which is determined by its i neighbouring particles within its smoothing 
length (hj). 

Using this concept, we can define the total mass (Mj) of a given k chemical element 
ejected by a j particle due to SNel or SNell, as 



where n v is the number of i neighbours of the j particle (n v ~ 40), mi and pi the gas 
mass and density of the i particle, respectively, Wiji — rj,hij) and hij, the symmetrized 
kernel and smoothing length, respectively. In this scheme, each i neighbouring particle gets 
a contribution of each metal defined by 



The j particle where the stellar mass has just formed is also included in this process. 



= ^m^MjWfo - rj , K,), 



(2) 



Ml = m i /p i M$W(r i -r j ,h ij ). 



(3) 
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The masses of H and He are proportionally decreased according to the metal mass received 
by each particle, so that its total mass is always conserved. 

The distribution of metals by using the SPH technique allows to spread out metals 
and enrich gas particles that may have not experienced SF but are close enough to other SF 
regions, working as an effective mixing mechanism. Note that gas particles can be enriched by 
more than one neighbour at each time-step. It has to be mentioned that numerical resolution 
would affect the mixing of metals since it also affects the determination of the neighbours. 
In Section 3, we will discuss this point in more detail. In the present implementation, metal 
mixing occurs only tided to the SF process, in the sense that, only new-born elements are 
distributed. 

One of the advantage of this implementation is that particles move according to the equa- 
tions of gravitation and hydrodynamics, leaving behind the closed-box hypothesis, among 
others. Particles with different astrophysical and chemical properties are mixed, particularly 
during violent events such as interaction and mergers, so that reproducing observations 
becomes a more challenging task. 

3 TESTS 

We carry out a series of simple tests to assess the effects of both, numerical resolution and the 
variation of the different parameters (C, 0sn, £sni) on the chemical properties of baryons. 
We also intend to evaluate the dependence of the results on the two adopted nucleosynthesis 
models, P98 and WW95. For this purpose, we follow the evolution of a homogeneous gaseous 
sphere. The sphere is represented by 10 3 gas particles initially distributed following a density 
profile p (x r~ 2 , in a radius of 625 kpc, without a dark matter halo. This simple experiment 
gives a clear idea of how each parameter affects the chemical properties of baryons and have 
the advantage of a low computational cost. 

3.1 Effects of numerical resolution 

Numerical problems may affect both, the process of star formation and specially the distri- 
bution of metals. With respect to the former, if dark matter haloes in galactic objects are 
numerically well resolved, the star formation process may not be strongly affected by the 
numerical resolution (see discussion in Dommguez-Tenreiro, Tissera and Saiz 1998). This is 
due to the fact that dark matter determines the potential well on to which the gas settles 
on. If dark matter haloes are well resolved, baryons are forced to distribute adequately. In 
this way, the gas density profiles are well reproduced and the SF can be correctly followed 
(Tissera 2000). Nevertheless, the gas component has to be resolved, at least, by few hundred 
particles (Navarro and White 1994). In our simulations, baryonic and dark matter compo- 
nents are represented by particles of equal mass. Therefore, dark matter haloes are resolved 
by a much larger number of particles than the baryonic component. Hence, the potential 
wells of the larger galactic haloes are very well defined ( see Tissera and Dommguez-Tenreiro 
1998). Consequently, we restrict our global analysis to galactic objects resolved with more 
than 250 baryonic particles. 

On the other hand, the radius of the neighbouring sphere in which the metals produced 
by the stellar component are distributed is affected by the numerical resolution of the experi- 
ment. Consequently, this problem may have a non-negligible effect on the metal distribution 
and the chemical properties of gas and stars, and requires a more detailed analysis. For 
this purpose, we have studied the evolution of two gaseous spheres without dark matter 
haloes. Each of them represents a system with total gaseous mass of 1.3 x 10 11 Mq and 
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initial radius of 615 kpc. One of the experiments (El) resolves the sphere with 10 3 particles 
(m p = 1.3 x 10 8 Mq), while the second experiment (E1H) has higher numerical resolu- 
tion: 8 x 10 3 gas particles (m p = 1.7 x 10 7 Mq). Since these spheres have no dark matter 
haloes, taking into account the former discussion, it is reasonable to expect that their star 
formation histories will be different. In fact, the sphere in E1H starts forming stars before 
its counterpart in El. However, the effects that this earlier SF has on their metallicity is 
not dramatic, as it can be observed from Fig.l, where we have plotted the mean values of 
[O/Fe] vs. [Fe/H] for both experiments. The greatest differences are obtained for the low 
metallicity stars ([Fe/H] < —3) due to the fact that star formation begins before in E1H 
and has a greater intensity than SF in El. As stars with higher metallicities are formed, 
this difference diminishes. Note that error bars in E1H comprise the distribution of El for 
[Fe/H] > —4 (except for the first point). 

The main difference observed between these two experiments is that low numerical res- 
olution produces a reduction in the dispersion of the distribution as can be seen from the 
error bars in Fig.l. This fact can be understood by considering that, in the experiment with 
lower resolution, the neighbouring sphere (n v ~ 40) comprises a larger volume, thus produc- 
ing a more effective mixing of the metals in the gas. Conversely, when numerical resolution 
is increased, the new elements produced by a certain particle are distributed in a smaller 
volume. This leads to a less efficient mix at the scale of the system, producing an increase in 
the metallicity dispersion. We cannot theoretically assess at what extent this dispersion has 
any physical meaning since the mixing mechanisms in galaxies are not yet well understood. 
It might be possible that SN explosions spread metals in larger volumes than that provided 
by the SPH technique. This fact would produce a more uniform-metal distributions in the 
interstellar medium, as it actually occurs in the low numerical resolution run. 

Taking into account these experiments, we conclude that the main effect of using the 
SPH technique to distribute metals is that the better resolved objects (i.e., those with higher 
number of baryonic particles that, in our simulations, implies larger masses) would have a 
larger dispersion in their chemical properties than the smaller ones. Mean values of the 
quantities will be considered as adequate estimations. 

3.2 Effects of the variation of the star formation parameters 

We are now going to analyse the effects that the different chemical model parameters have 
on the properties of the stellar populations of the spheres. 

As can be seen from Table 1, experiments El, E2 and E3 differ among each other in the 
life-time of binary systems that end up as SNel. The values considered are 10 8 , 4 x 10 8 and 
10 9 yr, respectively. The SF efficiency and ©sn used are the same in these three experiments. 
The ratio [O/Fe] vs. [Fe/H] clearly illustrates the effects of varying the parameters. Fig. 2 
shows [O/Fe] vs. [Fe/H] for these three experiments at different stages of evolution. As 
can be seen from this figure, as soon as the first stars are formed, the differences among 
them are very important. The smaller t^m used in El allows a rapid production of Fe and, 
consequently, the number of low metallicity stars ([Fe/H] < —3) is considerably lower than 
in E2 and E3, for all times. The larger tsNi value adopted in E3 produces a dex difference 
in the ratio [O/Fe] compared to the values in El and E2 for t = 1.13 x 10 10 yr. As time 
goes on, all experiments reach solar values, although, the variation of tsm among them leads 
to appreciably different evolutionary paths. For El and E2, the [O/Fe] is above observed 
values for the solar neighbourhood, and the steepness of the relation is much abrupt for 
[Fe/H]> —1. Experiment El (tsm = 10 8 yr) gives a more consistent relation for stars with 
[O/Fe] < 1. In this case, the final shape and range of abundances values are set as soon as 
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Table 1. Test: Main Parameter 
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WW95 


E2 


5e-6 
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WW95 


E3 


5e-6 
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10 


WW95 


El 


5e-6 


4 


1 


WW95 


E5 


5e-2 


2 


1 


WW95 


E6 


5e-6 


2 


1 


P98 



Units: [t SN i] = 10 8 yr ; [C] = Mpc 9 / 2 / Mg 2 /yr. 

the first stars formed which implies that even at early stages of evolution, this combination 
of SN parameters can produce stellar populations with metal contents up to solar. 

In Fig. 3, we show the same relations for experiments E4, E5 and E6. In E4 we have 
changed the relative ratio ©sn from 2 (in El) to 4. Since SNRII is the fraction of stars with 
masses greater than 8 M in a certain time-step, given by the IMF adopted, an increase 
in the ratio ©sn implies a decrease of SNRI. The effect of the variation of this parameter 
gives the expected differences in the abundance tracks. Note that there is no change in the 
slope of the [O/Fe] vs. [Fe/H] relation, but the abundances in E4 do not get to solar values 
because there are not enough SNel to produce the necessary amount of Fe. 

In experiment E5 we use a higher SF efficiency than that in El. In this case, stars are 
formed faster and the gas is enriched very efficiently so that the new-born stars immediately 
evolve to higher metallicity regions. As a consequence, the low-metallicity tail in the [O/Fe] 
vs. [Fe/H] relation vanishes. Note that the higher the C value used, the larger the stellar 
masses (A star ) formed, so that the number of SF events in a given baryonic particle dimin- 
ishes. This is the reason why the number of points in E5 is smaller when compared to the 
distribution in El. The important aspect to remark is that, in E5, the range of metallicity 
covered by the stellar populations is considerably smaller. This implies that a very efficient 
SF process can enrich the medium to solar values very quickly so that most of stars would 
tend to be located in the high metallicity region. 

The comparison between experiments El and E6 clearly illustrates how the simulation 
results depend on the two yield models, P98 and WW95. The features of the [O/Fe] vs. 
[Fe/H] relation in E6 can be appreciated in the bottom raw of Fig. 3. The main characteristic 
observed for this experiment, where P98 yields were used, is the lack of the low-metallicity 
tail ([Fe/H] < —3), which can be appreciated in El. On the other hand, the low-metallicity 
stars in E6 ([Fe/H] < —2), have almost constant values of [O/Fe], in contrast with the results 
of El, where a gradient is observed. This difference arises because P98 yields produce more 
iron with respect to the other elements than WW95, up to an order of magnitude more, in 
some cases. So, stars are rapidly brought to lower abundances ratios with respect to iron. 
Despite this fact, these relations are quite similar, with exception that WW95 seems to reach 
slightly higher metallicities values for the same model parameters at z = 0. This fact implies 
that P98 yields for enriched stars produce, in general, more elements than those of WW95. 
The comparison of other chemical elements shows that the major difference is in the [C/Fe] 
vs. [Fe/H] relation. In Fig. 4, we show this relation for experiments El and E6 as a function 
of time. As clearly seen in the case of E6, there is an excess of carbon in relation to iron for 
[Fe/H] > —1.2 that produces a positive gradient. 

To sum up, we find that all SF and SN parameters have non-negligible effects on the 
chemical properties of baryons. Then, in order to choose the correct combination of SN and 
SF parameters, a detailed comparison with observations has to be done. 
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4 GALAXY FORMATION 

We are now interested in applying this model to the study of galaxy formation in a cos- 
mological framework. As a first step, we look at the global properties of galaxy-like objects 
(GLOs) and compare them with observations. The different formation and evolutionary his- 
tories of each GLO (i.e., collapse time, merger tree, properties of progenitors, interactions, 
etc.) may affect their SF and chemical content in a complex way so that, even in the same 
experiments, GLOs may exhibit different characteristics. This fact makes of cosmological 
simulations a very useful tool for chemical evolution studies. 

The traditional model to study chemical evolution is the so-called One-Zone Simple 
Model based on two main assumptions: the system is isolated (i.e, no inflows or outflows) 
and well-mixed (i.e., instantaneous recycling), at all times (van der Bergh 1962; Schmidt 
1963; Tinsley 1980). Other hypothesis are related to the initial condition of the gas, IMF 
and nucleosynthesis yields. We remark the first two assumptions since they are inconsistent 
with the formation and evolution of galaxies in a hierarchical clustering scenario. In our 
experiments, structures form within a cosmological model suffering physical processes such 
as mergers, encounters, inflows, etc., which may affect the dynamics and kinematics of the 
dark and baryonic matters. Although some authors have included inflows and relaxed the 
hypothesis of instantaneous recycling (e.g., Chiappini et al. 1997), these models are not 
formulated in a cosmological framework and do not include dynamical effects that may 
affect the SF process and the mixing of chemical elements. Regarding the IMF and initial 
condition of the gas, we adopt a Salpeter IMF for all times, and the gas is assumed to be 
initially in primordial abundances. 

4.1 Numerical Experiments 

We performed SPH simulations consistent with a Cold Dark Matter (CDM) spectrum with 
Q, = 1, A = 0, Q h = 0.1, and a 8 = 0.67. We used N = 262144 particles (M part = 2.6 x 
10 8 M Q ) in a comoving box of length L = 5/i _1 Mpc (H = 100/i" 1 km s" 1 Mpc -1 , h = 0.5), 
starting at z — 11. Note that dark matter and baryonic particles have the same mass. The 
gravitational softening used in these simulations is 3 kpc, and the smaller smoothing length 
allowed is 1.5 kpc. Simulations SI to S5 share the same initial conditions while S6 shares 
the SF and SN model parameters of S2, but has different random phases. 

Simulations include SF and metallicity effects as described in Section 2. According to the 
discussion carried out in the previous Section, the value of C, Gsn and tsNi affect the chemical 
composition of the stellar and gaseous components. In order to assess the impact of these 
parameters on the chemical properties of galactic objects, we have performed simulations 
with the same initial condition but varying the model parameters (see Table 2). We also 
compare results from simulations where the two adopted yields, P98 and WW95, have been 
used. We will focus on the study of global chemical properties of the ISM and the stellar 
populations of GLOs, and their relation with the dynamical parameters of the objects. 

In order to study the properties of GLOs at z — 0, we identify them at their virial radius 
(Sp/p ~ 200; White and Frenk 1991). In Table 3 we give their total dark matter (Adark) and 
baryonic (Abar) number particles within the virial radius and their virial circular velocity 
(Kir)- Letter b in the label code of the GLO (second column, Table 3) indicates if the main 
baryonic system is formed by a pair of galactic objects. ^From the set of GLOs identified, 
we are only going to analyse those with more than 250 baryonic particles within their virial 
radius. 

As it is well known, GLOs are formed by a dark matter halo that generally hosts a main 
baryonic clump and a series of satellites. All quantities measured at the virial radius are 
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Table 2. Cosmological Simulations 



s c 


©SN 


*SNI 


yields 


SI 5e-6 


2 


1 


P98 


S2 5e-6 


2 


1 


WW95 


S3 5e-4 


2 


1 


WW95 


S4 5e-4 


3 


5 


WW95 


S5 5e-5 


2 


1 


WW95 


S6 5e-6 


2 


1 


WW95 


S7 5e-6 


2 


1 


WW95 


Units: [tsm] 


= 10 8 


yr; [C] 


= Mpc 9 / 2 /MQ /2 /yr. 



related to these complex systems. However, if we want to confront the simulated results 
with observations, it has to be taken into account that observed astrophysical quantities 
come from the luminous matter. For this purpose, we define a galactic object (GAL) as the 
structure determined by the main baryonic clump (including the dark matter mixed within 
it) hosted by a GLO. 

The radius that encloses 83% of the luminous mass of an exponential disc corresponds to 
the isophote of 25 mag arcsec -2 . Assuming the mass-to-luminosity ratio to be independent 
of radius, we define the optical radius of a GAL as the one that encloses 83% of its baryonic 
mass. All chemical and astrophysical properties will be referred to two optical radius (2R opt ), 
unless otherwise stated (see Table 3). This definition allows us to carry out a more meaningful 
comparison with observations. In this respect, V opt , defined as the circular velocity at 2R opt , 
is determined by the baryonic and dark matter distributions in GALs, and consequently, 
it is related to the physical mechanisms responsible for the concentration and distribution 
of the matter in the central region of a GLO (i.e., star formation, mergers, gas inflows, 
environment). Conversely, V v - W is a global parameter determined by the total potential well 
of the system. 

Recall that, in these simulations, there are three types of baryonic particles: pure gaseous, 
hybrid and total stellar ones. We identify those that belong to a given GLO and its GAL, 
and look at their chemical properties within 2i? opt at z = 0. 

4.2 Global Chemical Properties 

In this Section we study the correlations between chemical and dynamical properties of 
the GALs, in order to explore the possible physical mechanisms that may determine the 
metallicity of a galaxy, and to assess the dependence on model parameters. 

Observationally, defining the metallicity of a galaxy is a complex matter since it depends 
on the type and quality of available data, the element used as the estimator, etc. (e.g., 
Kunth and Ostlin 1999). Moreover, observations only give information on the metallicity 
of certain regions within a galaxy so that they could be considered good estimators of the 
global metallicity only if the ISM were efficiently mixed. 

Let us first review the main hypotheses of the one- zone Simple Model in order to under- 
stand the differences and advantages of our chemical scheme. The two basic assumptions of 
the Simple Model are: a) the system is isolated and b) it is well mixed at all times. It is also 
generally assumed instantaneous recycling. 

Concerning hypothesis a), our galactic objects form in consistency with a hierarchical 
clustering scenario via the aggregation of substructure. At z = a galactic object generally 
consists of a dark matter halo, a main baryonic clump and a series of satellites. Even if the 
objects do not have nearby similar-mass objects, they are rarely isolated. Moreover, there 
are always both, gas infall from the dark matter haloes and encounters with satellites that 
may tidally induce gas inflows as have been reported in numerical (e.g., Mihos & Hernquist 
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Table 3. Galaxy-like Objects: Main Parameters 
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Table 3 - continued 
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1996; Tissera 2000) and observational works (Barton, Geller & Kenyon 1999). At higher 
z, interactions and mergers notably increase so that the hypothesis of 'closed-box' is never 
valid since what is observationally identified as a galaxy is just a component of the whole 
system. It has also to be stressed that the star formation history is affected by the evolution 
of the galactic object. This evolutionary process determines the different chemical properties 
of the objects in our simulations (see Tissera et al. 2000 for details). 

Due to the mixing mechanism adopted which depends on the local gas density, and 
the fact that the systems are never isolated, the metallicity of the gas and the new-born 
stellar population is not uniform. At a given time, we found a significant spread in the 
values of the ISM metallicites within a given GAL (Tissera et al. 2000). These results are 
consistent with those from chemodynamical models (e.g., Samland et al. 1997). Regarding 
instantaneous recycling, we follow the evolution of the populations according to the stellar 
masses including supernovae I and II explosions. A consistent implementation of SNel is 
performed as described in Section 2.2. Planetary nebulae (PNe) have not been included in 
this work. 

Unfortunately we cannot treat in as much detail the ISM as chemodynamical models 
do, but this drawback is compensated by the fact that galaxy formation is well- described 
according to a cosmological framework without ad-hoc hypotheses. Moreover, the coupled 
non-linear evolution of dark matter and baryons has non-negligible effects on the star for- 
mation (e.g., Navarro & Steinmetz 2000; Tissera et al. 2000) that is directly related to the 
enrichment process. 

We will first attempt to assess which dynamical parameters correlate with the chemical 
properties of the GALs. We define global quantities for the chemical content of GALs. A 
global metallicity (Z) is assigned to the stellar population and to the gaseous component of 
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a GAL, taking into account all contributions from either, the stars or the gas respectively, 
according to 

where a refers to the component (i.e., stars or gas), Mk is the total mass of the k chemical 
element present in the a component, n is the total number of k chemical elements considered, 
and M a is the total mass of the a component within two optical radius. The estimated Z stax 
and Z gas for each GAL are listed on Table 3. 

We have analysed the GALs formed in simulations S2, S6 and S7. These experiments have 
the same star formation and cosmological model parameters but different random phases in 
the initial conditions. The SN parameters used in these experiments are taken from the best 
results given by the test runs (tsNi = 10 8 yr, 9sn = 2). In Fig.5a and Fig. 5c we show the 
global metallicities of the stellar population (Z star ) and the gaseous component (Z gas ) of the 
GALs versus their gas fractions. We find the expected trend indicating that the smaller the 
left-over gas, the higher the metallicity. Note that the slope of this relation for the gas and 
the stars is different, being steeper for the former. The averaged values of global metallicities 
for these runs are: < Z star /Z Q >= 0.25 ± 0.05 and < Z gas /Z Q >= 0.24 ± 0.09. Note that 
we are considering the metal content of stars and gas particles within 2i? opt irrespectively 
of their particular location. These mean metallicity values cannot be directly compared to 
that of disk stars in the solar neighbourhood. 

In Fig. 5b and 5d, we plot Z stSLT and Z gas versus the total stellar mass (M star ) within 
2R opt . As can be seen there is no clear trend between these parameters, albeit a weak 
tendency to have the higher metallicities in both, gaseous and stellar components of GALs 
with intermediate stellar masses: 3 x 10 10 M Q to 6 x 1O 1O M . 

For comparison, we have estimated these relations for our GALs assuming they have 
behaved according to the Simple Model. We include three cases: O.2Z (dotted-dashed lines), 
0.5Z Q (solid lines) and Z Q (dashed lines) Q . It can be appreciated the significant differences 
between the Simple Model and the results of the numerical simulations. The slope of Z stw jZ & 
versus M gas /Mb ar is different as well as the systematic increase of stellar metallicity as a 
function of the stellar mass formed in the Simple Model which contrasts with the behaviour of 
the numerical simulations. By inspection of the same relation for the gaseous components, we 
find that the Simple Model overestimates the metallicity of the galactic gaseous components 
for the parameters that match the chemical distribution of the stars. Hence, it seems that it 
is not possible to reproduce simultaneously the trends found in the stellar population and 
the gas of the simulated GALs with the Simple Model. This is somewhat expected since our 
code takes into account more complex processes (such as mergers, interactions, gas infall, 
etc) that are not considered in the Simple Model. 

We have studied the dependence of Z stav on V v - n and V opt finding no significant corre- 
lations. We expected the lack of correlation with V™ since it is determined by the whole 
gravitational bound system, while the chemical properties of the GALs are more strongly 
related to the fate of the baryonic matter in the internal regions. In fact, Fig. 6a shows no 
significant trend between the stellar mass of the GALs and their virial velocity. On the con- 
trary, there is a clear correlation of M star with V opt as can be seen from Fig.6b. If a universal 
value for the stellar mass-to-light ratio is assumed, this correlation would imply a Tully- 
Fisher relation similar to those found in previous works (e.g., Tissera et al. 1997; Navarro & 

* In Fig. 5 (b,d), lines represent the estimations of the Simple Model for a system with total mass Mt, ar = 10 10 Mq, initially 
with M s tar = and Z = 0. 
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Steinmetz 2000). Given that our GALs follow the Tully-Fisher relation, one would expect, 
in principle, a correlation of the stellar population metallicity with the V opt which is not 
found. We think that this lack of correlation is due to the fact that the chemical properties 
of these GALs are also affected by their merger history. Violent events have a strong impact 
on the baryon distributions and, consequently, on the metals mixing (White 1981; Cora et 
al. 2000). Furthermore, these objects are continuously accreating gas that contributes with 
pristine material to form new stars. From these results we find that the mean metallicity of 
the stellar population of the GALs cannot be directly linked to the optical mass either. Note 
the lack of SN wind effects in our simulations, as first discussed by Dekel & Silk (1986), could 
strongly affect the evolution of the gaseous component of low virial mass haloes, producing 
a non-negligible effect on the chemical properties of these systems. 

We have also studied how the properties discussed above change for different nucleosyn- 
thesis yield models. We have analysed the relations shown in Fig. 5 for experiment SI which 
has been run with the same SF and SN parameters than those used in S2, but with the 
nucleosynthesis model of P98 instead of those of WW95. It can be seen in Fig. 7a, a similar 
correlation than that shown in Fig. 5a although with a shallower slope; and a lack of corre- 
lation between Z staT and M star in Fig. 7b. It can also be appreciated that global metallicities 
in the simulation with P98 model are higher than those in simulations using WW95 yields, 
with average values of < Z gas /ZQ >= 0.40 ± 0.11 and < Z sta _ T /ZQ >= 0.45 ± 0.06. Hence, 
the different yields used imply significantly different results. 

In order to explore the effects of changing the model parameters, 0sn, ^sni and C, we 
plot in Fig. 8 Z stax versus M gas /M bar and M star for simulations S3, S4 and S5 (see Table 2). 
By inspection to Fig. 8a we can see that this correlation is present for all these experiments 
except S3. This simulation has the most efficient star formation, the largest value of C, which 
is an important parameter that strongly affect star formation and metallicity: for higher C 
values, and keeping the SN parameters fixed, the dispersion in metallicity increases. Note 
however that if the rate of SNel is decreased (larger 6sn) and the tsm increased, this 
correlation is recovered (S4). Hence, these three parameters are relevant at determining the 
global chemical properties of the GALs. In the case of the relation between metallicity and 
luminosity (total stellar mass), the correlation is not present as it can be seen from Fig.8b. 
However, we cannot further study the relative importance of the model parameters unless 
we look into more detail to the chemical properties of the gaseous and stellar components 
(to be analysed in Tissera et al. 2000). 

Estimates of the mean Z star and Z gas values of GALs in each simulation (see Table 4) 
show that when the gaseous component is more gradually transformed into stars (SI, S2, 
S6, S7), the mean global metallicity of stellar populations and the ISMs at z = are very 
similar, independent of the nucleosynthesis yield models used. Conversely, when the SF 
efficiency is increased (S3, S4, S5), the difference between the mean Z stlir and Z gas becomes 
very important. In this case, most of the metals are locked into stars, regardless of the SN 
parameters. Hence, in our chemical model a very efficient SF process produces, on average, 
ISMs considerably less metal-rich than the stellar populations of the GALs. 

However, it has to be noted that the mixing mechanism plays a key role in determining 
this result. A different implementation could lead to a more efficient distribution of metals 
in a shorter time-scale, quickly enriching the ISMs. Nevertheless, the way chemical elements 
are mixed in the ISM is still a controversial question that remains to be resolved from a 
theoretical and observational point of view. 

Unfortunately, we cannot directly compare these results with those derived from ana- 
lytical or chemodynamical model since the latter focus to the study of the Galaxy and do 
not have a sample where dependences with the dynamical parameters of the galactic objects 
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Table 4. Mean Z s tar and Z gas values of GALs 

SI S2 S3 S4 S5 S6 S7 

< Z gas /Z > 0.45 ±0.11 0.21 ±0.06 0.29 ± 0.14 0.18 ±0.06 0.21 ± 0.06 0.22 ± 0.06 0.27 ±0.10 
<Z sta . r /Z Q > 0.43 ±0.06 0.26 ±0.05 0.44 ± 0.15 0.25 ± 0.05 0.31 ± 0.06 0.24 ±0.04 0.26 ±0.06 

could be studied. Raiteri et al. (1996), although describing a similar chemical model as the 
one studied in this paper and using the SPH technique, do not consider different galactic 
haloes so that we cannot compare our findings with their results. 

4.3 Comparison with Observations 

Observations of galactic and extragalactic HII regions and OB associations provide infor- 
mation about the chemical properties of the forming stellar population and the ISM. These 
observations give relations for some primary elements such as S, Ne, and C, and for the 
so-called secondary ones, like nitrogen, as a function of the ratio O/H. In this Section, we 
resort to these observational results to assess the global chemical properties of the stellar 
populations and ISM of the simulated GALs. 

Concerning primary elements, Galactic and extragalactic HII regions show that the ratios 
(S/O) and (Ne/O) do not depend on (O/H). However, the dispersion in these ratios is quite 
large (Pagel 1997). The Simple Model actually predicts a similar behavior since in this 
scheme, the ratios between primary elements are constant. However, this model fails to 
reproduce the observed correlation between (C/O) vs. (O/H). A similar problem is detected 
for the observed (N/O) vs. (O/H) that exhibits a steeper correlation than that predicted by 
this model. Note that most observations are obtained from HII regions which are assumed 
good tracers of the metallicity of the galactic ISM. 

We analyse the ISMs in our GALs, which are determined by the gas particle properties 
within 2R opt - The global metallicity, Z^ as , for each k chemical element in the ISM of a GAL 
is defined as follows, 



1VJ.CT 



where n p is the total number of i gas particles within 2R opt , m\ is the mass of the k th element 
in the i th particle and M gas is the total gas mass within 2R opt . An equivalent relation can 
be defined for the stellar population, Z^ t3X . 

In Fig. 9, we plot log(S/0) for the ISM in the GALs in simulations S2, S3, S4 and S5. 
As it can be seen from this figure, the mean values obtained in the simulations are in very 
good agreement with observations (taken from Pagel 1997). The average log(S/0) for GALs 
shows a large dispersion that compares well with the observed values. The larger S/O ratios 
are obtained for the simulation with the highest SF efficiency (simulation S3). However, by 
changing SN parameters in order to diminish the effects of SNI explosions, and mantaining 
the same SF parameter in the simulation (simulation S4), we obtain substantially smaller 
abundance ratios and dispersions. 

In Fig. 10, we plot log(C/0) versus log(0/H) for GALs in simulations S2, S3, S4 and S5 
(WW95) and simulation SI (P98). We include observations of HII regions (Pagel 1997). A 
remarkable fact that can be seen from this figure is the large difference between the results 
of WW95 and P98 models. GALs in SI, which are exactly equivalent to those in S2, except 
for the nucleosynthesis model adopted, have ratios almost half an order of magnitude larger 
than their WW95 counterparts, and are out of the observational range. Conversely, GALs 
in any of the WW95 runs have ISMs with abundances similar to observed ones. However, 
they do not show the gradient reported from observations of Galactic and extragalactic HII 
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regions (e.g., Garnett et al. 1995; Kennicutt and Garnett 1996; Kobulnicky and Skillman 
1996), although the simulated abundances extend only in the range —3.8 <log(0/H) < —3.2. 
The simulated mean values are lower than the observed ones but it has to be recalled that 
we are not including PNe yields that are thought to be important contributors of C (and 
N). 

Observations of (N/O) in HII regions show that this ratio increases with (O/H) (e.g., 
Pagel 1992). A gradient in the secondary element is actually predicted by the Simple Model, 
although the steepness of the predicted relation is much larger than that actually observed 
in HII regions. There are several hypothesis that may explain the origin of these differences 
(Pagel 1997). We estimate this relation for GALs in WW95 runs. As shown in Fig. 11, the 
simulated ISMs have abundances that are in good agreement with observations. However, 
we only find GALs with 12+ log(0/H) < 8.5, while the observed gradient in HII regions 
is 12+ log(0/H) > 8.5. We do not find GALs with average ISM metallicity corresponding 
to solar abundances, neither extremely low metallicity IZW18 type objects. It should be 
considered, however, that only massive GALs are analysed in the simulations. 

Another important observational correlation is that between 12+ log(0/H) and the gas 
fraction estimated for irregular and blue compact galaxies (e.g., Axon et al. 1998). This 
relation can be used to assess the metal content of the ISMs in relation to the left-over gas 
mass of the corresponding GALs. In Fig.12, we plot 12+ log(0/H) versus log (M bar /M gas ) 
for the ISM in the simulated GALs (M bar and M gas are the total baryonic and gaseous mass 
within 2_R opt , respectively). Observations taken from Pagel (1997) have been included for 
comparison. The estimated abundance ratios in the models are consistent with observations. 
It is notable that the observed relation is well reproduced in the sense that objects with a 
lower gas fraction tend to have larger metallicities without the need to introduce other 
physical mechanisms such as SN energy injection effects to the ISM. The Simple Model 
predicts a decrease of the metallicity with the gas richness of the system but it fails to 
correctly reproduce the observational relation. 

5 CONCLUSIONS 

We have implemented a chemical model in a SPH cosmological code. First results are re- 
ported together with the assessment of the effects of model parameters and numerical reso- 
lution. It is found that the major problem introduced by a low numerical resolution is the 
artificial smoothing of the chemical properties of the objects. However if a minimum number 
of baryonic particles is imposed, then average values can be considered reliable. 

The mean metallicity is found nearly independent of the total stellar mass of the system 
indicating that the effects of mergers, interactions and gas infall on the mass distribution 
and metal mixing are very significant. We find correlations between the stellar mass M star 
and the virial and optical circular velocities (Kir, V opt ). However, the correlation with V opt 
is significantly tighter, consistent with the Tully-Fisher relation. Nevertheless, the global 
metallicity of GALs, Z star shows no dependence on V opt , in spite of the fact that these 
GALs satisfy the Tully-Fisher relation and shows a strong correlation between the global 
metallicity of their stars and gas components with the left-over gas fraction of the systems. 
We only find a significant correlation between the global GALs metallicity of both stars and 
gas with the left-over gas fraction of the systems. 

An analysis of the GALs assuming they have behaved according to the one-zone Simple 
Model show important differences with the results of the numerical simulations. The slope of 
Z star /Z versus M gas /Mb ar and the systematic increase of stellar metallicity as a function of 
the stellar mass formed in the Simple Model contrasts with the behaviour of the numerical 
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simulations. We find that the Simple Model overestimates the metallicity of the galactic 
gaseous components if it is constrained to match the stellar chemical content. According to 
the Simple Model, the global metallicity of a system continuously increases with the stellar 
mass in disagreement with the results found for the simulated GALs that show higher 
metallicities for intermediate mass objects. 

We find that the star formation efficiency, SN model parameters and the nucleosynthesis 
yields significantly affect the chemical properties of the GALs. In those models where the 
SF process is gradual the mean stellar and gaseous metallicities are similar. 

The observed abundance ratios for primary and secondary elements in Galactic and 
extragalactic HII regions are naturally obtained. However, the ranges of the simulated global 
metallicities are smaller than the observed ones. Also, the relation between the (O/H) and gas 
fraction obtained in the numerical simulations is consistent with the observational results. 

The suitable agreement between the models and the observations for these relations 
suggests that hierarchical clustering scenarios are able to reproduce the chemical properties 
of galaxies. This is a very encouraging fact taken into account that galactic objects formed by 
the accretion and mergers of substructures. Violent events are common and ubiquitous, and 
can affect the dynamical evolution of the matter, regulating the SF and chemical evolution 
in galactic objects. 

We intend to improve this chemical model in the future by allowing the gaseous com- 
ponent to cool according to its metallicity, and a further study of the mixing process of 
chemical element in the ISM. It is also under study a more efficient decoupling mechanism 
between stars and gas that could shorten the period of hybrid state of a baryonic particle. 
Finally, energy feedback remains an open question that we hope to address in a future work. 
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FIGURE CAPTIONS 

Figure 1: Comparison of the [O/Fe] versus [Fe/H] relation obtained from the low (solid line) 
and high resolution tests (dashed lines). 

Figure 2: [O/Fe] vs. [Fe/H] for experiments El, E2 and E3 where the life-time of binary 
systems has been varied: 10 8 yr, 4 x 10 8 yr and 10 9 yr, respectively. Time t is given in units 
of 10 10 yr. 

Figure 3: [O/Fe] vs. [Fe/H] for experiments E4, E5 and E6 where the relative ratio of 
SNII/SNI (Osn), the star formation efficiency (C) and the nucleosynthesis yields have been 
changed with respect to El. Time t is given in units of 10 10 yr. 

Figure 4: Comparison of the [C/Fe] vs. [Fe/H] relation for experiments El: WW95 and 
E6: P98, yields. Time t is given in units of 10 10 yr. 

Figure 5: Global metallicities of the stellar population and gas component in GALs as a 
function of their gas fraction, M gas /M bar , (a,c), and their total stellar mass ,M star , (b,d) in 
simulations S2 (open pentagons), S6 (filled squares) and S7 (open squares). Lines represent 
the relations given by the Simple Model for 0.2Z Q (dotted-dashed lines), 0.5Z & (solid lines) 
and Z Q (dashed lines). M star is given in units of 10 10 M Q . 

Figure 6: The total stellar mass of GALs (M star ) as a function of a) the virial velocity 
(Kir) and b) the optical velocity (K pt ) for the same GALs shown in Fig. 5. Velocities are 
given in kms -1 and M star is given in units of 10 10 M Q . 

Figure 7: Global metallicities of the stellar populations in GALs (Z staT /Z Q ) in simulation 
SI (P98) as a function of a) the gas fraction (M gas /M bar ), and b) the total stellar mass 
(M star ) of GALs. M star is given in units of 10 10 M . Lines represent the relations given by 
the Simple Model as in Fig.5. 

Figure 8: Global metallicities of the stellar populations (Z star /Z Q ) as a function of a) the 
gas fraction (M gas /Mb ar ), and b) the total stellar mass (M star ) of GALs in simulations S3 
(filled pentagons), S4 (open triangles) and S5 (open stars). M star is given in units of 10 10 M Q . 

Figure 9: (S/O) versus (O/H) for the ISM of GALs in simulations S2, S3, S4 and S5 (see 
Fig.5 and 8 for feature code). Observations of HII regions have been included (small filled 
circles). 

Figure 10: (C/O) versus (O/H) for the ISM of GALs in simulations S2, S3, S4 and S5 (see 
Fig.5 and 8 for feature code). GALs in SI (filled triangles) have been included. Observations 
of HII regions have been included (small filled circles). 

Figure 11: (N/O) versus (O/H) for the ISM of the GALs shown in Fig. 9. Observations 
of HII regions have been included (small filled circles) . 

Figure 12: Logarithm of the oxygen abundances versus the logarithm of their gas fraction 
for the same GALs shown in Fig.9. Observations of HII regions in blue compact and irregular 
galaxies from Pagel (1997) have been included (small filled circles). 
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